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1  Introduction 


The  preparation  of  an  a  priori  prescribed  finite-dimensional  unitary  transformation  via 
the  evolution  of  an  externally  driven  quantum  system  is  a  theoretical  and  technological 
challenge  which  plays  an  important  role  in  several  fields  including  atomic  and  molecular 
manipulation  [1,  2,  3,  4],  quantum  computation  [5,  6],  quantum  cryptography  [7].  This 
paper  will  describe  a  systematic  technique  to  generate  arbitrary  special  unitary  transfor¬ 
mations  in  a  quantum  system  which  draws  on  concepts  from  control  theory.  The  results 
described  here  are  applicable  to  the  broader  domain  of  quantum  system  manipulation. 

The  principal  concept  exploited  here  is  that  any  special  unitary  matrix  can  be  decom¬ 
posed  into  products  of  transformations  with  a  particular  structure.  Similar  decompositions 
have  recently  been  suggested  as  well.  Reck  and  Zeilinger  [8]  use  the  decomposition  of  any 
m  xm  unitary  matrix  into  a  product  of  simpler  ones  (i.e.,  tensor  products  of  an  arbitrary 
2x2  block  and  a  complementary  block  consisting  of  the  identity  matrix  of  dimension 
(m  -  2)  obtained  in  [9]  to  produce  an  optical  implementation  of  quantum  cryptography 
schemes.  In  any  factoring  approach  one  first  needs  to  determine  the  appropriate  structure 
within  the  factors  and  then  show  that  any  special  unitary  matrix  can  be  written  as  a 
product  of  such  factors.  In  this  regard,  Eberly  and  Law,  [10],  considered  a  decomposition 
for  controlling  the  quantum  state  of  a  cavity  field.  There  is  a  relationship  between  their 
factors  and  those  in  the  present  paper,  and  the  number  of  factors  needed  are  different. 
Finally,  DiVincenzo  and  collaborators,  have  demonstrated  that  any  unitary  matrix  can  be 
written  as  a  product  of  unitary  matrices  which  are  either  a  tensor  product  of  an  arbitrary 
2x2  block  and  Jm-2  or  a  special  4x4  block  and  Im-\  [11].  This  latter  decomposition 
was  motivated  by  a  multiparticle  implementation  of  quantum  logic  gates.  Here  we  will 
consider  the  fundamental  case  that  the  non-trivial  sub-block  is  of  size  2x2.  Unlike  the 
conventional  decomposition  where  the  2  x  2  blocks  are  in  U(2),  the  approach  presented 
in  this  paper  respects  the  selection  rules  appropriate  to  the  system.  Moreover,  it  can  be 
tailored  to  take  into  account  practical  laboratory  constraints. 
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The  specific  ingredients  of  our  approach  to  creating  the  unitary  matrices  has  the  fol¬ 
lowing  steps: 

•  Employ  an  atom  or  other  simple  quantum  system  whose  energy  levels  are  well  sepa¬ 
rated,  such  that  there  are  no  interfering  resonances  with  any  pair  of  optically  coupled 
levels.  Assurance  of  this  condition  may  call  for  the  introduction  of  suitable  static  ex¬ 
ternal  fields.  Spectral  separation  permits  the  system  to  be  controlled  by  sequentially 
addressing  pairs  of  levels.  Furthermore,  this  circumstance  allows  one  to  employ  the 
rotating  wave  approximation  when  considering  control  field  designs.  Controllability 
is  also  possible  under  conditions  more  relaxed  than  this  clean  spectral  separation, 
but  each  such  case  must  be  individually  analyzed  for  its  applicability. 

•  Decompose  the  prescribed  NxN  unitary  matrix  into  a  product  ofiVxiV  unitary 
matrices  which  are  non-trivial  only  in  a  sub-block.  Here  we  will  consider  the  sim¬ 
plest  case  where  the  sub-block  is  of  size  2x2.  This  decomposition  will  respect  the 
selection  rules  appropriate  to  the  system,  and  will  also  be  tailored  to  take  into  ac¬ 
count  the  rotating  wave  approximation  and  any  limitations  on  the  matrix  elements 
of  the  electric  dipole  operator,  which  is  assumed  to  be  the  coupling  operator.  The 
decomposition  used  in  this  paper  is  different  from  the  standard  one  where  the  2x2 
block  can  be  any  element  of  U( 2)  (see  [9]). 

•  An  explicit  characterization  of  the  exponential  of  an  arbitrary  real  linear  combination 
of  the  Pauli  matrices  will  be  utilized.  This  formulation  will  be  useful  in  parametrizing 
the  2x2  blocks  of  the  previous  steps  to  provide  an  explicit  design  for  the  laboratory 
control. 

The  balance  of  this  paper  is  organized  as  follows.  Section  2  presents  a  derivation  ofthe 
decomposition  needed  for  the  physical  system.  Theorem  2.1  demonstrates  that  any  matrix 
in  SU( 2)  can  be  decomposed  as  desired.  In  addition,  it  is  shown  how  selection  rules  can 
be  respected  during  the  decomposition  process.  Remark  2.1  explains  how  the  number 
of  factors  may  be  minimized.  Setion  3  discusses  how  the  concepts  in  the  paper  may  be 
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applied  to  state  preparation  and  observation.  Finally  the  Conclusion  section  4  addresses 
the  extent  to  which  the  ideas  in  this  paper  may  be  generalized  to  other  applications. 

2  Generation  of  Atomic  Scale  Unitary  Transforma- 
tions 

Consider  an  atom  irradiated  by  a  seuqence  of  electromagnetic  pulses  with  the  aim  of 
guiding  its  evolution  initially  from  the  identity  matrix  to  a  desired  final  unitary  matrix. 

The  atom  is  assumed  to  have  the  following  properties: 

1.  The  atom  has  M  accessible  levels  which  are  well  separated  for  controllability  based 
on  frequency  discrimination. 

2.  The  matrix  elements  of  the  dipole  operator  (with  respect  to  the  basis  consisting  of 
the  eigenfunctions  of  the  free  Hamiltonian)  are  all  real  and  the  diagonal  elements 
are  identically  zero.  This  latter  assumption  can  be  relaxed  in  certain  cases. 

3.  The  atom-pulse  pair  is  “controllable”  in  the  following  sense.  There  is  always  a 
connection  between  any  two  levels,  i.e.,  if  selection  rules  preclude  accessing  level  i 
directly  from  j,  then  there  is  a  set  of  levels  (possibly  more  than  one)  forming  a  ladder 
such  that  state  i  can  be  accessed  at  least  indirectly  from  state  j. 

The  multi-level  atom  will  be  irradiated  by  a  sequence  of  tailored  pulses.  Each  of  the 
pulses  above  will  address  only  a  prescribed  pair  of  levels.  There  is  flexibility  in  the  pulse 
structure  and  each  should  be  chosen  to  be  as  simple  as  possible.  The  action  of  the  entire 
sequence  of  pulses  is  described  below. 

The  evolution  of  the  unitary  generator  is  described  by  : 

iU{i)Q)  =  (H0-  nc)U(t,Q)  (2.1) 

where  is  the  induced  dipole  operator  and  e(t)  is  the  control  field.  Defining  Cl(t,  0)  by 

U(t,  0)  =  exp(~iHot)Q(t,  0) 
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leads  to  the  interaction  representation 

iO(<,  0)  =  -e(t)eiHo,ne-iHotn(t,  0) 

The  control  pulse  is  taken  to  be  e(<)  =  A(t)cos(ujt  +  <j>).  where  <j>  is  a  phase  to  be  chosen 
and  A(t)  is  assumed  to  be  slowly  time  varying  compared  to  w_1,  and  the  frequency  u  is 
resonant  with  levels  i  and  j  to  be  coupled. 

The  M  x  M  matrix  €(t)eiHoi  pe~iHot  will  have  all  allowed  couplings  given  by  but, 
given  the  conditions  stated  earlier,  we  may  neglect  all  terms  exccept  the  (i,j) th  and 
(J,  i)th  entries  which  respectively  are  pijel*A(t)  and  mje~t<f>A(t)  (recall  that  all  entries 
of  the  matrix  representation  of  /i  are  assumed  to  be  real  and  the  diagonal  entries  are  all 
zero). 

The  pulse  is  applied  for  a  time  tx  with  the  only  restriction  that  the  integral  f*1  A(t)dt 
be  insensitive  to  a  slight  change  in  tx.  Applying  the  rotating  wave  approximation  yields: 

n(ti,0)  =  71,0(0,0))  =  IdM 

where  the  matrix  V\  is  a  M  x  M  matrix  that  is  the  identity  except  for  a  2  x  2  block,  which 
is  the  exponential  of  a  2  x  2  matrix  of  the  form 


7i  is  the  complex  conjugate  of  the  complex  number  71  which  has  the  following  polar 
representation: 

71  =  Hije'*  f  A(t)dt  (2.3) 

Jo 

It  will  be  seen  later  (  Equation  (2.5))  that  the  magnitude  of  the  complex  number  71 
enters  into  our  considerations  only  as  the  argument  of  trigonometric  functions.  Hence 
A{t)  and  <j)  are  chosen  so  that  7i’s  magnitude  can  take  any  desired  value  within  [0,2tt]  to 
achieve  the  necessary  action  to  generate  Vi.  This  point  is  useful  for  practical  laboratory 
considerations. 

After  the  application  of  this  first  pulse  the  unitary  generator  is  given  by: 

U(h,0)  =  e-iHot'V 1 
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Now  another  pulse  of  length  t2  is  applied  which  is  resonant  with  only  some  other  allowed 
pair  of  levels.  We  then  have: 

U(t,0)  =  U(t,h)U(h,0) 

Once  again  setting  (Mi)  =  gives: 

iQ(Mi)  =  -e(t)eiHo^~tl)ne~iHoii~u)Q(t,ti) 

Premultiplying  both  sides  of  the  last  equation  by  eiHoU ,  using  the  rotating  wave  approxi¬ 
mation  and  then  integrating  up  to  time  t  =  t 2  yields: 

eiH°ua(t2,tl)  =  V2eiH°un{tl,tl) 

Here  V2  is  again  a  M  x  M  matrix  which  is  the  identity  except  for  a  2  x  2  block  which  is 
the  exponential  of  a  2  x  2  matrix  similar  in  structure  to  that  in  (2.2).  As  ft(<Mi)  is  the 
identity  matrix,  the  last  equation  becomes: 

fi(fa.*i)  =  e-iHotlV2eiHoh 

Hence,  U(t2M)  =  e~iH^-u^{t2M)  =  e~iH^V2eiHat\  and 
U(t2, 0)  =  £f(*s,<i)£f(*i,0)  = 

If  a  total  of  k  such  coherently  locked  pulses  are  applied  (starting  at  t  =  0  and  lasting 
for  a  total  period  of  time  T  =  <i  +  t2  +  . . .  +  tjt)  we  get: 

U{T,  0)  =  e-iHoTVkVk-l...V2V1  (2.4) 

where  the  Viyi  =  are  M  x  M  matrices  with  each  being  the  identity  except  for 

a  2  x  2  block  which  is  the  exponential  of  a  2  x  2  matrix  similar  in  structure  to  that  in 
(2.2),  with  7i  replaced  by  suitable  complex  numbers  7 j.  If  it  is  desired  that  U(T,  0)  be 
a  prescribed  matrix,  V,  then  the  question  that  needs  to  be  addressed  is  whether  V  can 
be  decomposed  as  in  the  last  equation  by  a  suitable  choice  of  pulses.  Equivalently,  since 
e-*H0T  -ls  a|ready  known,  the  question  is  whether  any  M  x  M  unitary  matrix  can  be 
decomposed  as  the  product  nf_1  V{  for  some  choice  of  k?  If  the  answer  is  affirmative  then 
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one  can  consider  optimizing  the  sequence.  This  optimization  would  seek  to  minimize  k, 
while  simultaneously  considering  the  pulse  shaping  capabilities  of  the  apparatus. 

It  is  necessary  to  show  that  7  can  be  written  as  e^n^T*)-  Equivalently  we 
have  to  show  that  every  unitary  matrix  V  =  etHoTV  can  be  written  as  U^=lV(jk)-  It  will 
be  shown  below  that  this  can  always  be  achieved  if  V  is  in  SU(M).  As  any  W  G  U(M) 
can  be  written  as  eirV  with  T  a  real  scalar  and  V  G  SU(M),  this  shows  that  the  goal  of 
creating  V,  up  to  a  scalar  phase  factor  which  is  irrelevant,  has  been  attained. 

In  the  remainder  of  the  paper  we  will  assume  that  V  €  SU{M)  and  denote  by  VH  its 
Hermitian  conjugate.  Towards  the  end  of  demonstrating  that  V  has  the  desired  factoriza¬ 
tion,  we  premultiply  VH  by  a  product  of  matrices  with  each  factor  in  the  product  being  a 
matrix  which  is  the  identity  except  for  a  2  x  2  block  of  the  form  V(jk),  so  that  the  result 
is  the  identity  M  x  M  matrix. 

Before  proceeding  to  the  general  situation  we  first  illustrate,  via  examples,  how  selection 
rules  are  respected  in  the  process  of  reducing  VH  to  the  identity  by  left-multiplication. 
Assume  that  the  atom  has  M  =  4  levels  and  that  the  only  allowed  transitions  are  1  —  2 
and  1  — ►  3  and  1  -*  4.  Any  two  levels  are  still  accessible  from  one  another,  albeit  via  level 

1. 

Denoting  VH  =  (a,j),t,j  =  we  will  first  reduce  column  4  to  (0, 0, 0, 1)T, 

then  column  3  to  (0,0, 1,0)T,  then  column  2  to  (0, 1,0, 0)T  and  then  finally  column  1  to 
(1,0,0,0)T.  The  reduction  will  be  achieved  by  premultiplying  VH  by  suitable  matrices 
which  are  tensor  products  of  I2  with  a  V(yk)-  Notice  that  this  order  is  not  the  only  possible 
one.  All  that  is  important  is  that  column  1  be  reduced  last.  In  order  to  appreciate  this 
point  suppose  to  the  contrary  that  column  1  has  been  transformed  first  to  (1,0, 0,0)  , 
thereby  leaving  a  matrix  of  the  form 


1 

0 

0 

0 

0 

ai 

6i 

Cl 

0 

t>2 

C2 

l 0 

&3 

C3 
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with  the  3  x  3  matrix  on  the  bottom  right  corner  an  element  of  5(7(3).  Once  this  is 
achieved  we  will  not  be  able  to  reduce  the  remaining  columns  (if  any)  to  unit  columns 
because  level  1  will  no  longer  have  any  population  in  it  to  effect  any  further  intermediate 
transitions  required  between  the  remaining  levels. 

Returning  now  to  the  procedure  suggested  above,  consider  the  reduction  of  column  4 
to  the  corresponding  unit  column.  We  may  view  the  fourth  column  as  a  particular  state 
of  a  4  level  system.  First  transfer  all  of  the  amplitude  from  the  a2 4  element  to  the  ai4 
element  by  using  the  1  -+  2  transition.  This  can  be  effected  by  premutiplying  VH  by  a 
matrix  in  5(7(4)  of  the  form: 

/  \ 

a\  &i  0  0 

Ci  d\  0  0 

0  0  10 

^  0  0  0  1  j 


The  submatrix 


bi 


l  is  in  5(7 ( 2)  whose  elements  encode  information  about  024 


y  cx  di  j 

and  aH  (see  Remark  2.1  below).  In  the  process  of  this  premultiplication,  VH  will  be 
transformed  into  a  special  unitary  matrix  with  its  (2,4)  entry  equal  to  0.  Next  use  the 
1  -+  3  transition  to  transfer  all  population  from  the  (3,4)  entry  of  the  new  matrix  to  its 
(1,4)  entry.  This  maybe  effected  by  premultiplying  the  new  matrix  by  a  matrix  in  5(7(4) 
of  the  form: 

02  0  62 
0  1  0 
C2  0  d-x 
^000 


\ 

0 

0 

0 


The  submatrix 


(  a2 


in  5(7 ( 2)  encodes  information  about  the  (3,4)  and  (1,4) 


y  °2  fa  j 

entries  of  the  new  matrix,  as  before.  This  premultiplication  will  result  in  a  5(7(4)  matrix 
whose  (2,4)  and  (3,4)  entries  are  both  zero.  Finally  the  1  — ►  4  transition  is  used  to 
transform  the  the  fourth  column  of  the  last  matrix  to  the  column  (0,  0, 0, 1)T.  Once  again 
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this  may  be  achieved  by  premutiplying  the  last  matrix  by  a  matrix  in  SU{ 4)  of  the  form 


az  0  0  63 
0  10  0 
0  0  10 
C3  0  0  dz 


l 


where  the  matrix 


0E3  63 

C3  dz 


in  SU(T)  is  specifically  tailored  to  achieve  the  desired  effect. 

After  this  sequence  of  premultiplications  the  net  effect  is  a  matrix  in  SU{ 4)  whose  last 
column  is  (0,0,0, 1)T.  Unitarity  requires  that  the  4th  row  be  (0,0,0, 1).  Special  unitanty 
now  forces  the  remaining  3  x  3  block,  V,  to  be  in  SU( 3).  We  now  proceed  to  reduce  the 
third  column  of  V  to  (0,0, 1)T  by  making  use  of  transition  1  —  2  to  transform  all  the 
amplitude  to  level  1,  and  then  using  the  transition  1  —  3  to  reduce  the  third  column 
of  the  resultant  matrix  to  (0,0, 1)T.  As  before  both  of  these  transitions  are  achieved  by 
premultiplication  by  matrices  in  517(4)  whose  last  column  and  last  row  are  (0,0,0, 1) 
and  (0,0,0. 1).  The  balance  of  these  matrices  are  tensor  products  of  I\  and  the  matrix  in 
517(2)  specifically  tailored  to  achieve  the  desired  effects. 

Remark  2.1  Let  us  consider  the  structure  of  the  matrices  in  517(2)  which  achieve  the 

desired  effects  above.  For  example  in  nulling  out  the  024  element  of  we  have  several 

a14 

choices  for  the  matrix  in  517(2).  Let  r  be  the  norm  of  the  vector  ,  and  then  seek 

(a24  / 

a  matrix,  U12  €  517(2)  such  that 

V12 

Usually  9  can  be  taken  as  arbitrary.  Only  in  the  last  step  in  converting  any  column  to  the 
corresponding  unit  vector  does  0  have  to  be  zero.  Fixing  9  leads  to  a  unique  V12  which 
will  achieve  the  above  result.  This  is  just  the  statement  that  517(2)  acts  transitively  on 
spheres  in  C2  of  any  given  radius.  However,  the  freedom  in  the  choice  of  0  (except  in  the 
last  step)  is  useful  in  minimizing  the  number  of  V(jkYs  needed.  Indeed  it  will  be  shown 


\ 

it 

\ 

ai4 

re'6 

024  J 

° 

/ 
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below  that  any  matrix  in  SU( 2)  is  a  product  of  at  most  3  matrices  of  the  form  V{lk  )■  So 


is  a  matrix  of  the  form  V(y)  such  that  ^(7) 


\ 

/ 

V 

— 

\ q } 

!  \ 

v'IpP  +  MV"  ^ 


0 


where 


/ 


the  complex  numbers  p  and  q  have  polar  representations  p  —  |  p  |  e  **  and  q  \  q  \  e  ^ . 

( p  ^ 

Notice  that  the  phase  of  p  is  the  same  as  that  of  V(y) 


.  This  will  be  demonstrated 


w 

in  Lemma  (2.1)  later.  Thus  not  insisting  that  9  be  zero  (except  in  one  step)  reduces  the 
number  ofV(yk)’s  to  used  and  thus  the  number  of  pulses  to  be  applied. 


The  general  case  of  an  M-level  atom  can  be  handled  by  directly  extending  the  above 
procedure  to  express  any  matrix  in  SU(M )  as  a  product  of  matrices  which  are  tensor 
products  of  IM-2  and  an  element  of  SU( 2).  Thus  all  that  remains  is  to  show  that  every 
element  of  517(2)  can  be  written  as  a  product  of  V(7*)’s  and  this  is  achieved  by  the 
theorem  below. 


Theorem  2.1  Any  matrix  V  €  SU( 2)  can  be  expressed  as  the  product  nj^Vfr*)-  Fur- 
thermore,  l  can  be  taken  to  be  no  more  than  three . 


Proof:  The  proof  proceeds  via  the  following  steps: 


1.  Suppose  7  =  rel<i>  is  given.  Then 


W(j)  =  exp  ( 


0 


17  \ 


0  J 

Let  R?  =  a2  +  /J2  (thus  R  =  2 r).  Employing  the  formula  for  the  exponential  of  an 
arbitrary  real  linear  combination  of  the  Pauli  matrices  given  in  the  Appendix,  with 
6  =  0,  7  =  we  get: 


V(y)  = 

^  ™  R 
COS  2 

\ 

0 

2  sin  4 

4-  1 

^  R 

/ 

0 

a-ifi  ^ 
2 

0 

cos  f  ) 

2 

<x+iP 
^  — 2~ 

0  y 

(  R 
cos  f 

\ 

0 

zsin  r 

0 

\ 

re'* 

0 

COS  J  J 

+ 

r 

re-^ 

°  j 
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Hence, 


V(y)  = 


cos  r  ie1^  sin  r 


j'c-,^sin  r  cos  r 


(2.5) 


2.  For  5  e  St/ (2),  the  following  holds: 

/ 


S  = 


e,Acosa  etA*sino 
e*(7T— /i)  sin  Q,  e”,Acosc* 


sir 


with  A  >  0  and  ft  >  § ,  and  some  real  a 

l  .  \ 

a  b 


Proof:  If  V"  = 


c  d 


,  then  the  rows  of  of  K  are  orthogonal  and  detV  —  L  So, 


if  a  =  cos  ae,A  and  6  -  sinae*^,  then  7  admits  the  above  description.  Here  A  and 
fi  are  chosen  to  satisfy  the  inequality  constraints  required.  Clearly  A  can  always  be 
chosen  positive.  If  does  not  already  satisfy  the  required  inequality  constraints, 
then  we  may  write  b  =  e*(7r+*i)  sin(— a),  and  then  take  tt  4-  p  as  the  new  /z,  and  —a 
as  the  new  a.  Since  the  original  /i  could  always  be  chosen  positive,  this  achieves  the 
desired  inequality  constraint  on  p.  Also,  as  cos(— a)  =  cos(a),  a  is  unaffected  by 
this  process. 

3.  The  number  of  required  for  writing  5  as  a  product  of  them  can  be  taken  to 

be  no  more  than  three. 


Proof:  Let  71  =  aet9 ,  where  6  —  A  +  /j  -  Then  S  =  W(ji) 


l  , 

e 


i\  0 


\ 


V 


0  e 


-*A 


/ 


^  eiA  0  ^ 


Now 


v 


=  W{ 72)1^(73)  where  72  =  fe"71  and  73  =  \tx°2  where  07 

0  e-A 

and  02  are  any  real  numbers  such  that  A  =  o\  —  <Ji  +  Jr.  This  completes  the  proof  of 


Theorem  2.1. 


Finally,  the  lemma  which  was  used  in  Remark  (2.1)  follows: 
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\ 

P 

w 


in  C 2  there  exists  a  complex  number  7  such  that: 


Lemma  2.1  Given  a  vector 


W{  7) 

Proof:  Choose  7  =  aeie ,  where  cos  a  =  and  9  =  p,-  A  -  §. 

Clearly 


\ 

\ 

p 

,  - _ 

n/I  p  I2  + 

q  |2e^ 

K  9 ) 

V 

/ 

Thus, 


sin  a 


9  I 


IpI2  +  M2 


w{rt) 


\ 

v 

9  / 


/  \ 

■  \ 

cos  a  isinae*^  A  2' 

\p\e'» 

i  isin  2 )  cos  a  j 

{  |9|€‘A  ) 

^  Ippe^  +  ilgpe'^-f)  ^ 


0 


/ 


VT  P  I2  +  U  |V'M 
0 


....  \ 

/ 


3  Applications  to  State  Preparation  and  Observation 

This  section  will  consider  applications  of  the  generation  of  special  unitary  matrices  to  a) 
state  preparation,  and  b)  state  observation.  Case  a)  will  be  discussed  from  the  perspective 
of  controllability,  [12]  and  case  b)  aims  to  show  that  the  preparation  of  certain  sequences 
of  special  unitary  matrices  leads  to  the  determination  of  the  state  of  the  system. 
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•  State  Preparation:  The  case  of  state  preparation  is  implicitly  present  in  a  previous 
paper  [12],  where  the  problem  of  controllability  of  molecular  systems  was  analyzed 
by  considering  the  system  defined  by  the  corresponding  unitary  generator.  However, 
a  slight  refinement  is  needed  to  utilize  that  analysis,  since  in  this  paper  we  have  onlj 
shown  how  special  unitary  matrices  may  be  prepared.  Given  two  vectors  u  and  v 
on  the  unit  sphere  of  CM  there  is  an  entire  family  of  unitary  matrices  U  such  that 
Uu  =  v.  To  construct  them  we  first  consider  the  unitary  matrices  Uu  and  Uv,  where 
Uu  and  Uv  have  as  their  first  columns  the  vectors  u  and  v.  The  remaining  columns 
can  be  any  vectors  which  render  Uu  and  Uv  unitary  (e.g.,  those  obtained  by  applying 
the  Gram-Schmidt  process).  Clearly  Uup  =  u  and  Uvp  =  v  where  p  =  (1, 0, . . . ,  0)r. 
Then  the  matrix  U  =  UVU?  satisfies  Uu  =  v. 

Now  consider  the  case  where  the  matrix  acting  on  u  to  give  v  is  desired  to  be  special 
unitary.  Construct  Su  =  AU “  and  Sv  =  UVB  according  to  A  =  diag(l,  e‘s\  1, . . . ,  1), 
where  ei6‘  =  detUu,  and  B  =  diag(l,e-<4*,  1, . . 1)  with  e1'4’  =  det Uv.  Then  both 
Su  and  Sv  are  special  unitary  and  thus  the  product  5  =  SVSU  is  also  special  unitary 
and  satisfies  Su  =  v,  since  Uu=v.  Hence  by  preparing  the  special  unitary  matrix  S 
and  acting  on  the  initial  state  u  it  is  possible  to  create  the  desired  state  v.  In  most 
situations  u  =  (1,0, ..  .,0)T. 

It  would  be  interesting  to  examine  this  approach  for  initial  state  creation  for  low 
dimensional  systems.  The  ability  to  create  initial  states  for  low  dimensional  systems 
can  be  combined  with  other  approaches  such  as  optimal  control  [2,  13,  3,  4]  to  create 
desired  states  for  higher  dimensional  systems. 

•  State  Observation:  State  observation  is  fundamental  to  control  engineering  [14], 
and  it  is  equally  significant  for  quantum  systems  [15].  The  ability  to  determine 
that  a  desired  state  has  been  actually  obtained  is  needed  to  assess  the  success  of 
the  control  process  and  its  use  in  a  variety  of  applications.  We  assume  that  the 
measurement  process  can  determine  the  population  of  some  particluar  level  N.  Thus, 
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if  ’/'Final  =  EZi  cie<</>i&’  where  i  =  ■  •  •  -  M  is  some  basis  (e'S’  eigenfunctions 

of  the  free  Hamiltonian),  then  we  are  able  to  measure  |  cN  |2.  By  preparing  the 
special  unitary  matrix  which  results  in  the  exchange  of  populations  between  levels 
N  and  n,  it  is  possible  to  measure  |  cn  |2,  for  all  n  =  1, . . M.  The  only  remaining 
quantities  to  be  determined  are  the  phases  <f>i,i  =  Since  only  relative 

phases  are  relevant,  we  assume  that  the  reference  phase  <j>N  is  0.  Next  the  special 
unitary  matrix  which  is  the  identity  except  for  the  following  2x2  block  in  the  (N,  n) 
position  is  prepared: 

1  1 
-1  1 

When  applied  to  Spinal  results  in  a  vector  which  is  identical  to  ^Final  except 
that  cN  and  c„  are  replaced  by  (,v  =  c-^  and  C»  =  s^-  By  assumption,  we 
can  readout  |  Cat  |2-  Since, 


V2 


V 


\ 

/ 


I  C N  P  =  ^{1  CN  I2  +  |  Cn  I2  +CN Cn  +  C*NCn} 

......  •  ,,  ICw|3-7(|cn13+|c»|2)  _  - 

it  is  possible  to  determine  Mn  —  - 2|caH  |c»l  ”  2  ” 

cos (<t>N  -  <f>n)-  However,  this  does  not  determine  4>n  since  the  cosine  function  is 
double  valued  on  (0,27r). 

To  remedy  this  latter  problem  we  prepare  the  special  unitary  matrix  which  is  the 

tensor  product  of  Im— 2  and  the  following  2x2  block  in  the  (IV,  n)  position. 

/  \ 


y/2 


1  i 
i  1 


When  applied  to  Spinal  ’  tbe  transformation  results  in  a  vector  which  is  identical  to 
Spinal  except  that  cn  and  c„  are  replaced  by  Cn  —  CN^a‘  and  Cn  =  •  Once 

again,  by  assumption  we  are  able  to  measure  |  Cn  P-  Since  the  quantity 

Mn  =  - - - =  Sin (<pN  -  <t>n) 

can  be  measured,  it  is  clear  that  {<j>N  -<t>n)  can  be  determined  for  each  n,  and  hence 
so  can  <j>n,n  —  1, . . .,  M. 
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A  total  of  3(Af-l)  measurements  are  required  on  a  sequence  of  controlled  systems  for 
ihe  2(Af-  1)  unknowns  |  c,-  | 2,  The  above  scheme  is  a  special  case  of  determining 
the  elements  of  a  density  matrix  via  projection  measurements  developed  in  [15]. 

4  Conclusions 

In  this  paper  a  constructive  scheme  for  the  creation  of  any  special  unitary  matrix  has  been 
analyzed.  Although  the  physical  system  used  here  is  a  single  particle  interacting  with  an 
electromagnetic  field,  there  are  some  general  features  to  the  scheme  which  are  relevant  for 
other  quantum  systems  [16].  These  features  are: 

♦  The  quantum  system  in  question  can  be  controlled  by  addressing  subsystems  of 
dimension  two  sequentially. 

•  The  logarithm  of  the  unitary  generator,  obtained  after  the  controlling  mechanism  is 
applied  to  any  of  these  two  dimensional  subsystems,  lies  in  the  span  of  any  two  of 
the  Pauli  matrices. 

The  first  point  above  suggests  that  we  represent  our  desired  special  unitary  matrix  as  the 
products  of  factors  which  are  non-trivial  only  in  a  2  x  2  block.  If  no  special  structure  on 
this  2x2  block  is  desired,  then  this  can  always  be  achieved,  [9].  If  a  special  structure  is 
desired,  then  a  preliminary  question  is  what  special  structures  are  feasible?  The  second 
point  yields  the  answer.  The  reason  is  that  the  Lie  algebra  generated  by  any  two  of  the 
Pauli  matrices  equals  the  entire  Lie  algebra  of  SU( 2).  Thus,  if  the  second  point  is  valid 
for  the  quantum  system  being  studied  then  any  matrix  in  SU(2)  can  be  decomposed,  in 
principle  at  least,  into  a  product  of  matrices  with  the  prescribed  special  structure.  The  very 
important  practical  questions  of  obtaining  this  decomposition  in  an  algorithmic  fashion 
and  also  minimizing  the  number  of  factors  in  the  decomposition  depend  on  the  particular 
system.  However,  if  the  quantum  system  is  such  that  in  the  second  point  above  the 
logarithm  belongs  to  the  span  of  ic rx  and  iay,  then  the  constructive  algorithm  introduced 
here  is  applicable  to  it  as  well. 
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Although  quantum  computation,  [5,  6],  was  not  the  primary  focus  of  this  work,  its 
needs  provide  useful  background.  Any  proposal  to  fabricate  a  quantum  computer  encom¬ 
passes  three  steps:  i)  Preparation  of  the  initial  state  which  will  serve  as  the  input,  ii) 
Preparation  of  the  unitary  transformation  which  will  function  as  the  logic  gate;  and  iii) 
Reading  the  final  output.  With  every  concrete  protocol  proposed  for  a  quantum  computer 
the  physical  means  for  achieving  these  3  steps  must  be  specified.  In  addition  the  nature 
of  each  step  depends  on  the  particular  problem  or  quantum  algorithm.  This  paper  consid¬ 
ered  a  single  atom  (i.e.,  a  particle  in  the  limit  of  one  electron)  as  a  realization.  There  are 
many  practical  limitations  associated  with  this  choice  (as  is  the  case  with  any  other  choice 
thus  far  presented),  but  there  are  also  strong  arguments  to  suggest  that  an  atom  would 
be  an  ideal  system  for  practical  laboratory  study.  Single  particles  are  attractive  from  the 
point  of  view  of  isolating  the  system  from  external  influences,  and  thereby  enhancing  the 
lifetime  (provided  states  of  sufficient  radiative  lifetime  are  involved).  Furthermore,  both 
the  theoretical  and  experimental  aspects  of  controlling  molecular  and  atomic  dynamics 
are  advancing  very  rapidly.  Ail  three  steps  above  can  be  subsumed  into  the  problem  of 
creating  a  prescribed  unitary  matrix  from  the  evolution  of  a  controlled  atom.  It  is  pos¬ 
sible  to  achieve  entanglement  within  a  single  molecule  by  simultaneously  manipulating 
two  degrees  of  freedom  such  as  rotation  and  vibration,  or  radial  and  angular  degrees  of 
freedom  in  a  single  atom  etc.,  A  further  critical  matter  concerns  the  structure  of  the  algo¬ 
rithms  used  in  quantum  computation.  The  current  algorithms  are  designed  to  work  with 
multiple  particles.  However,  this  circumstance  does  not  preclude  working  with  multilevel 
systems  The  nature  of  the  state  observation  is  also  very  important.  Section  3  considered 
full  state  observation,  but  other  less  demanding  observations  are  desirable  in  keeping  with 
the  quantum  computation  algorithms  proposed  so  far.  These  topics  go  beyond  the  scope 
of  this  paper,  but  the  work  here  should  provide  motivation  for  exploring  these  questions 
further. 

Finally,  it  would  be  interesting  to  examine  the  field  of  molecular  control,  [2,  13]  j  from 
the  perspective  of  this  paper.  Most  problems  in  molecular  control  can  be  cast  as  state 
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preparation  for  mutilevel  systems.  Section  3  showed  that  this  goal  can  be  subsumed  into 
the  problem  of  creating  one  out  of  a  family  of  special  unitary  matrices.  The  main  feature 
of  our  algorithm  is  that  it  does  not  require  any  costly  computations,  although  the  pulse 
sequenxce  still  needs  to  be  determined. 


5  Appendix:  Exponential  of  Pauli  Matrices 

Since  no  reference  was  found  for  the  result  below,  a  proof  is  being  included  for  the  sake 
of  completeness.  It  facilitates  obtaining  decompositions  of  unitary  matrices  adapted  to 
various  constraints  in  the  control  pulse. 


Theorem  5.1  Let  ax  =  ~  \  I ,  =  \  I 

y i  o  y  \i  o 

the  Pauli  matrices.  Then  we  have 


and  crz  =  | 


0 


^P)/2+_A=sin(^I  )(a(t)<rx+0(t)<Ty++m<r> 


expi(a(t)<rx+p{t)<ry+6(t)ff*)  =  c°s(Xy"1)^2+ sin(— Y~ 

where  a(t),f3(t)  and  8(t)  are  real  valued  functions  of  time ,  and 


\(t)  =  (a(t)7+m2+m2) 


Proof:  First,  a  direct  calculation  shows  that  cr^  =  <Ty  =  -  \h\  Vx<?y  - 

—  Gy  (Tx ;  ffy&z  =  —  —GzCy]  and  that  CT  Z(J  x  —  —  —  CTX(TZ. 

Hence, 


(ia<rx  +  ifoy  +  iSazf  =  (ia<rx  +  i/3<ry  +  i6<rz)(iaax  +  i@<ry  +  i6tr. 


■[a2crl  +  ctf3ax<Ty  +  a6ax<rz  +  a0<ry<rx  +  /?2<r2  4-  /36trycrz  +  otb<rz<rx  +  60crzcry  +  b2a: 

oo  i  _»  00  i  _\ 

=  E  (2^r(— /2)" +  £  (to+T)f(  Th)n{iaa* +  if*'  + 

n=0  v  ' 


„e0(2n  +  D!v  4 


or  +  /r  +r)h 


Hence 


ei(o/.+/»/,+7/.)  =  g  J^y(ialx  +  ifil,  +  iyl z )2"  +  g  J^TTy{iaIx  +  ifiIy  +  ilh 
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-  E  + E  iSTipiiW”1- + W + -V. ) 

=  E  ^<T>2"'’  +  iS(toTI)!(T)i"(o/' + "• +  rt  > 


1  ,  -A 


n=0 


=  + th  +  7/. ) 


=  cos(\/A)/2  +  ^^sin(^~ )(a/x  4-  &Iy  4* 7^)- 


which  is  the  desired  result. 
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